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Abstract 

The rate constants for recombination and exchange processes are studied in 
terms of two different flux correlation approaches: one is the Yamamoto approach, 
which is based on the hnear response theory, and the other is the Miller one. Using 
those approaches we consider two exactly solvable cases, i.e., the free particle and the 
parabolic potential models. Since the rate constants for recombination and exchange 
processes are calculated by Laplace transforms of the flux correlation functions, the 
two approaches give different results. In the present calculation, we find that the 
rate constant in the Yamamoto approach is larger than that in the Miller approach 
by about 40% at low temperature (~ 100 K) and high pressure (~ 1 GPa). The 
difference is considerable in the region where quantum effects dominate. 
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When one wants to study a chemical reaction at the most detailed level, it is neces- 
sary to calculate the Schrodinger equation for a state-to-state differential scattering cross 
section, which is a function of total energy E and total angular momentum J. Such 
quantum reactive scattering cross sections have actually been studied for simple chemical 
reactions, where a time-dependent scattering formalism based on the ^-matrix Kohn vari- 
ational approach^ or a coupled channel method in hyperspherical coodinatesi has usually 
been used. 

However, in chemical applications, there are many cases where only the (microcanoni- 
cal or canonical) rate constant for a reaction is needed. If the full, state-to-state scattering 
calculation has been carried out, the rate constant is, of course, given by an average of 
the cross sections. If it is, however, only the rate constant that is desired, such a complete 
calculation for all state-to-state information is not economical. Furthermore, because of 
the rapid growth of the number of open vibration-rotation channels with increasing ther- 
mally accessible collision energies, the calculation of a rate constant via exact quantum 
state-to-state calculations would not be feasible even for a simple reaction. 

The traditional way of evaluating a rate constant is the transition state theory (TST).i 
However, TST is a classical theory and it is approximate, because it does not involve the 
effect of recrossing of the system over the transition state dividing a potential surface. A 
number of improvements have been proposed to take the recrossing effect into account .0 i 

In the early 60's, Yamamotoi first formulated an exact expression for the rate constant 
as an application of the general statistical mechanical theory of irreversible process, which 
was established by Kubo et al.i and Mori.i Later (in the early 70's), Miller et al.li i 
separately developed a method for the rate constant using a time integral of the flux- 
flux autocorrelation function, which is exact in the limit that the dynamics is extended 
to t oo. Because the fiux-fiux correlation is calculated via time- dependent quantum 
mechanics, the feasibility of this approach depends on how to evaluate the Hamiltonian, 
flux and time evolution operators for the system. The fiux-fiux autocorrelation function 
method has been applied to a variety of chemical reactions.0 In particular, the reaction 
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oi H + H2 has been studied intensively.lill 

It has recently been shown how a quantum mechanical version of the Lindemann 
mechanism for coUisional recombination 

A + B ^ AB*, (1) 
AB* + M AB + M, (2) 

can be handled by the flux-flux autocorrelation function for the A — B collision.li3 
Here the process is affected by the bath gas M. Some applications of this new theory are 
listed in Ref.l4. It is furthermore possible to generalize the formalism to include chemical 
reactions as well as recombination: 

A + BC ^ ABC* -^AB + C, (3) 
ABC* + M ABC + M. (4) 

Equations (H) and simultaneously describe the recombination process {A + BC 
ABC) and the exchange reaction {A + BC AB + C)M This method has been applied 
to the interesting (combustion) reactions (O + OH ^ H + O2) and the recombination 
reactions {0 + OH + M —>■ HO2 + M <— H + O2 + M). Those reactions are very important 
in atmospheric chemistry.S 

The purpose of this paper is to study the difference between the way proposed by 
Yamamoto,i which is based on the linear response theory or the so-called Kubo formula,^ 
and the fiux-fiux autocorrelation function method proposed by Miller et alJ'i The two 
approaches can provide the same result to the rate constant for a simple chemical reaction, 
because it is given in terms of the integral of the fiux-fiux correlation function with respect 
to time. However, the shapes of the correlation functions calculated by the two methods 
are quite different from each other. For the recombination and exchange reactions (like 
eqs.(|]) and (|D), the rate constants in the two approaches could be different because they 
are evaluated by Laplace transforms of the fiux-fiux correlation functions. It is expected 
that the difference will appear in the region where quantum effects dominate. 
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In this paper, we first review the correlation function method briefly and show the 
difference between Yamamoto's and Miller's approaches explicitly. The rate constants for 
recombination and exchange processes are also discussed. Next, we study two exactly 
solvable cases, i.e., the free particle and the parabolic potential models. Finally, the 
summary and conclusion are given. 



1 Flux-flux correlation approach to rate constants 

In the classical limit, a rate constant is generally given by an average of the flux through 
some dividing surface that separates reactants from products (see Fig.|l|). The canonical 
rate constant is then given byEi 

KiiT) = Qr{T)-\2Ti)-f j dpi J dqi e-^^(^^^'^^^)F(pi,gi)P(pi,g-l), (5) 

where = ksT (T, temperature) and provides the initial conditions of the 

momenta and (reaction) coordinates for classical trajectories of the system (consisting 
of / degrees of freedom). The system is described by the Hamiltonian H{pi,qi). Here 
Qr{T) is the partition function per unit volume for the noninteracting reactants and F is 
the flux factor which describes the trajectories crossing the dividing surface specified by 
s{q) = 0: 

F{p,q) = f^h{s{q)) = S{smv,, (6) 

where s{q) is some function of position q that is negative on the reactant side and positive 
on the product side. Then, h{s) is the step function, which is +1(0) for s > (<)0, and 
Vg is the normal component of the velocity to the dividing surface s{q). The factor V in 
eq.(^) involves all information of the dynamics and it is unity when the trajectory is on 
the product side in the infinite future and zero otherwise. This implies that it is given by 

V{pu gl) = lim h{sm)) = / dt -/i(s(g(t))) = / dt F{p{t), m), (7) 

Jo dt Jo 

^We use the natural unit, i.e., h/2Tr = c = 1. 
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where eq.(||) is used. Thus, V provides the probabihty that the trajectory hes on the 
product side of the dividing surface at t — >^ oo. The rate constant then reads 



roo 

QrKl{T) = / dt Cclit) 

Jo 



where 

C4t) = (27r)-^ J dpi J dqi e-'5^(^^i''^i)F(pi, q,)F{p{t), (9) 

This means that the rate constant is calculated by the time integral of the flux-flux 
autocorrelation function Cdit). 

To take quantum effects into account, it is necessary to replace the phase space integral 
by a quantum trace representation. In the linear response theory,0 the response functon 
is usually defined as 

0BA(t) = -^tr(p[A,5(t)]), (10) 

where p = e~^^ /tY{e~^^) is the density operator for an equilibrium state. After the 
perturbation by the operator A at t = 0, the response of the quantity B{t){= e'^^Be"'*^) 
at time t in the system is described by the response function (pBAif)- For the fiux-fiux 
autocorrelation function, one can identify that A = h{s) and B{t) = F{t). Here F{0) is 
the fiux operator at t = 0, which is given bj 



F{0) = t[H,his)] = ^ 



m m 



with p the momentum operator and m the reduced mass of the system. Thus, the response 
function for the rate constant is 

<l)Fh{t) = -ttTip[h{s),Fit)]). (12) 

Using the Kubo identityB 

[A, e-^"] = e-^^ r dX e^^ [H, A]e"^^, (13) 
Jo 



the response function reads 



<pFh{t) = t d\ tT{pF{-iX)F{t)), (14) 







where eq.(|n]) is used and F{-i\) = e^^F{0)e~^^. 

In the hnear response theory, the relaxation function ^ba is defined as 

roo 

$sA(t)=hm/ ds(PBA{s)e-''. (15) 
The relaxation function for the rate constant is thus given by 

roo roo rf3 

<^Fh{t) = / ds (f)Fh{s) = ds d\ tT{pF{-tX)F{s)), (16) 
Jt Jt Jo 

where we assumed that the response function decreases rapidly as t — >■ oo. The rate 
constant in quantum mechanics is now given in terms of the relaxation function at t = 

roo rf3 roo 

Qrk{T) = p-^<l>FhiO) = / dt d\ tr(pF(-a)F(t)) = / dt C{t), (17) 

JO Jo Jo 

where the flux-flux autocorrelation function in quantum mechanics C{t) is deflned by 

C{t) = f dX tr{e-^^ F{-iX)F{t)). (18) 
Jo 

In the present notation, C{t) in eq.(|l^) corresponds to the flux correlation proposed by 
Yamamoto.i We should note that there exists an integral with respect to A which stems 
from the Kubo identity and that it is dispensable in the classical limit /3 — > 0. Because 
it is more convenient to use the commutation relation in eq.(|T2|) rather than eq.(|TB|) in 
actual calculations, we re-deflne the Yamamoto's correlation function by 

C^{t) = ^tT{e-^^[h{s),Fm, (19) 

where the superscript Y stands for "Yamamoto" . 

By contrast, in Miller's approach!' § the variable A in the flux is flxed to be P/2 and 
the A integral is performed. Thus, from eq.([l8|) Miller's correlation function is given by 

C^'^it) = tr(F(0)e***^F(0)e-'*^^), (20) 

where = t — i(3 /2 and the superscript M stands for "Miller". This modification certainly 
makes actual calculations simple, because the fiux operators are involved symmetrically in 
the correlation function. In fact, Yamamoto's correlation function C^(t) is not identical 
to C^^{t). However, their integrals with respect to time are identical to each other, and 
hence they can provide the same rate constant. Therefore, the way of Miller et al.li i 
certainly has some distinct advantages in actual numerical calculations. 



2 Rate constants for recombination and exchange 
reactions 



It is possible to generalize the flux correlation approach to treat recombination and ex- 
change reactions.0' 111 111 It may be intuitive and useful to begin with the classical de- 
scription of the process again. Let us consider the reaction of A + BC —>■ AB + C and 
ABC (see Fig.|^). The classical rate constants for the exchange {A + BC AB + C) and 
recombination {A + BC ABC) reactions are again given by eq.@, i.e., averages of the 
flux and the probability V{pi,qi) over the Boltzmann distribution. Here is 

the flux at the reactant dividing surface Sr (see Fig.^: 

Fr = —h{Sr) = 5{Sr)Vr. (21) 
(JjL 

Note that h{sr) is again the step function, which is 0(1) for position q to the left (right) of 
the dividing surface s^, and that Vr is the normal component of the velocity to the surface 
Sr- Similarly we deflne the step function for the product dividing surface Sp by h{sp) (see 
Fig.^]), that is, h{sp) = 0(1) for position q to the left (right) of the dividing surface Sp. 
The difference of those step functions, hc{q) = h{sr{q)) — h{sp{q)), is unity for position q 
between the two dividing surfaces (i.e., in the "compound" region) and zero outside. 

Because the probability of the system experiencing a deactivating {ABC* + M — 
ABC + M) collision with the bath gas M can be evaluated by 1 — e"''* at time t [r] 
describes the frequency of deactivating collisions and it depends on pressure P and T of 
the bath gas), the recombination probability is estimated as 

where r is the time the trajectory (it is on Sr aX t = 0) is in the compound region. Thus, 
using he and an integration by parts, Vrec readJll 

/•CXD rj /"OO 

Vrec= dthMt))-rA^-e-'^')= dt{e-^'-l){Fr{t)-Fp{t)), (23) 



10 dt JO 

where Fi{t) = hi{q{t)) (z = r or p) 



For the exchange reaction, the probabihty is given by e ^'^'^ , where Tp is the time the 
trajectory exists through the surface Sp. The probabihty is eventually obtained asEl 

roc fj roo 

Ve.c= rft(l-/ip(g(t))-e-^*= / dte-^'Fp{t). (24) 
JO at Jo 

Inserting those probability functions into eq.(|^), the rate constants for the recombina- 
tion and exchange reactions are given by 

QrKnT.P) = I dte'^^\CfXt)-Cfp{t)), (25) 

/■oo 

QrkTiT.P) = I dte-^'C:l{t), (26) 

where 

C^Ut) = (27r)-^|t/pi|rfgle-^^(^^-^^^)F,,(pl,gl)F,,(p(t),g(t)), (27) 
Cfpit) = {2T,)-f j dp, j dq,e-P''^^^''l^'^Fr{p,Ai)Fpm.m)- (28) 

Here the relation 

POO roo 

dt CfM = I dt C:lit) (29) 

holds because in the limit ?7 — the recombination rate should vanish. 

The transcription of the rate constants to quantum mechanics simply involves replacing 
the classical correlation functions by their quantum mechanical counterparts. As in the 
classical case, the rate constants for recombination and exchange processes in quantum 



mechanics are thus given by 

POO 

Qrk^eciT^P) = I dte-^\Cjf{t)-Cjf{t)), (30) 

/•oo 

QrkTiT,P) = I dte-^'Cjfit). (31) 

Then, the fiux-fiux autocorrelation functions are given by 

Cl{T) = ltr(e-^^[/^(..),F,(t)]), (32) 
Crp{T) = ltr(e-^^ [/.(..), F,(t)]), (33) 
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in Yamamoto's approach, while 

C^!.{T) = tr(F,,(0)e***^F,(0)e-'*^^), (34) 
C,^(T) = tr(F,(0)e^**^F,(0)e-^*^^), (35) 

in Miller's approach. Note that Fi{0) (i = r or p) is the quantum mechanical flux, which 
is again given by eq.(|lTl) with 6{si) and h{si), instead of 6{s) and h{s). 

Because the rate constant is calculated by the Laplace transform of the flux-flux au- 
tocorrelation function, it is clear that the two approaches give different results. It is 
expected that they will coincide with each other in the classical limit, but the difference 
becomes large in the region where the A integration in eq.(|18D cannot be ignored. 

3 Numerical calculations 

In this section we calculate the (canonical) rate constants for recombination and exchange 
reactions using the Feynman path integral technique.0 A huge calculation is usually re- 
quired to obtain the exact matrix elemens of propagators for a realistic system. Fur- 
thermore, it is necessary to consider some approximations and numerical techniques like 
Monte Carlo samplings^ to perform it. Because the aim of this paper is to show how the 
rate constant in the Miller approach is different from that in the Yamamoto case, it would 
be more intuitive and useful to consider a simple system rather than a complicated case. 
We here study two analytically solvable cases: i.e., the free particle and the parabolic 
potential models in one- dimension, and leave more elaborate calculations for nontrivial 
cases for a forthcoming paper. 

3.1 Free particle case 

We first study the free particle case (see Fig.^l). The propagator for the free particle in 
a coordinate representation can be easily calculated by the path integral .0 The matrix 
element of the flux operator Fi{i = r or p) in coordinate space is also found easily for the 
free particle system. For details, see Appendix A. 
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The flux-flux autocorrelation function in the Miller approach is eventually given by 



rp 



47r(t2 + /32/4)3/2 



(3 2mt^(P 
2^ + /?2/4 



exp 



2(t2 + /?2/4) 



(36) 



with d the distance between Sr and Sp (see FigJ^). Note that the correlation depends on 
only the distance d and is independent of positions Sr and Sp, as it should be. From this 
expression the correlation function C^'^{t) is easily obtained as 



(37) 



^7r(t2+/?2/4)3/2- 

Those correlation functions are illustrated by the dotted curves in Figs.^ and ^, in which 
we define Ci = md'^ /2f3 and take Ci = 1.0 to illustrate the correlation functions clearly. 
The rate constant for the reaction without recombination (i.e., in the limit rj ^ 0) can be 
obtained by the integral of eq.(^) with respect to time (see eq.(|lj 

1 



QrHT) 



ivrp 



(3^ 



By contrast, in the Yamamoto approach the flux- flux autocorrelation function is given 
by (see Appendix A) 



cut) 



1 



X 



27r/32^2t(t2 + /32) 



t2 + /?2 + X+{dt^ + (3^- tf' COS X 



exp 



mfSd"^ 



2(t2 + /52) 



^3/2, 



(39) 



where X = m/5^(i^/2t(t^ + jS'^). This is not identical to eg. (|36D . In particular, at short 
time it is divergent like ~ 1/ \^ although it is integrable. Note that it again depends on 
only the distance d. If we set Sr = Sp (or = 0), we obtain0 



cut) 



i^WTW-t) 



3/2 



(40) 



2V2t7T(3\t^ + f3^y/^' 

and, as expected, we can find that for the usual rate constant the Yamamoto correlation 
function gives Qrk{T) = l/27r/5, which is the same as that in the Miller approach (see 
eg. (PBD ) . Those correlation functions are shown by the solid curves in Figs.^ and|^. 

It can be seen from Figs.^ and |^ that the interference effect in the correlation is taken 
into account correctly in the Yamamoto approach (although the vibrating behavior is 



10 



inconvenient for numerical calculations) . Contrastingly, in the Miller correlation function 
the interference is averaged and the shape is quite smooth. Thus, it is very convenient for 
actual computation. For the usual rate constant, the two approaches certainly give the 
same result, as we have seen above. 

Next, we calculate the rate constants for recombination and exchange reactions. The 
rate constants are given by eqs.(^) and (|3ll) . If we define the (rr)- and (rp)-Ta.te constants 
by 



Qr^rr 



Y,M 
rr 



Qr ^r 



rp 



dt e-'^'C^;'\t), 



(41) 
(42) 



Y,M 
rec 



the rate constants are given as k. 
Then, the Miller approach gives 

1 



l,Y,M _ i.Y,M J /,y,M _ r.Y,M 
'^rr '^rp ""^^^ '^exc ~ '^rp ■ 



M 
rr 



M 
rp 



lTi(3 Jo 

1 



dx 



+ 1/4)3/2^ 



^7lp Jo 



dx e~ 



+ 



exp 



a;2 + 1/4 



(43) 
(44) 



(X2 + 1/4)3/2 (x2 + l/4)5/2_ 

where x{= t/j3) is a dimensionless variable and a = f3r], while in the Yamamoto approach 
we find 

L-x)3/2 



Y 
rp 



2s/2iT/3 Jo 

2V27r/3 Jo 



dx e~ 



dx 



x^ 



x(a;2 + 1) 



a;(a;2 + 1 
exp 



(45) 



X2 + 1 



X 



(Vx2 + 1 + x)^/2 sin X' + (Vx2 + 1 - x)3/2 X' 



(46) 



with X' = Ci/x{x^ + 1). 

In order to convert the collision frequency t] to more familiar variables, we approxi- 
mate the collisional deactivation rate constant by an expression given by the hard sphere 
collision theory. Furthermore, if one uses the ideal gas expansion, the frequency can be 
expressed by0 

7/ = k,,act[M] = Pi/^ X 10-'\ (47) 



11 



with f] in fs ^, P in Pa and T in K. Then, we find 



a = Pri- 3.24 x ^ x 10"^ (48) 

The factor Ci is also converted as 

Ci = = 0.0103 X ATd\ (49) 

with A the reduced mass of the system in atomic mass units and d in A. In this paper we 
consider a system which has a small reduced mass (like H + CO ^ HCO or H + O2 
HO2) to illustrate the difference between the two approaches clearly. In the following 
calculations, we thus take A = 2 and = 2A and vary T and P. 

Now we are in a position to show our results for the free particle case. First we define 
ratios 

Rrr{T,P)=kJjk^,, (50) 

Rr,{T,P) = kl/k^^. (51) 

Figures ^ and ^ illustrate the two ratios for the free particle case. Here we choose T = 
100 ~ 400 K and P = 0.1 ~ 1 GPa. (To check the accuracy of the present numerical 
calculation, we have also evaluated the rate constant for the free particle with rj = 
and compared the result with the exact value given by eg. (|55|) . It is confirmed that the 
numerical calculation is sufficiently accurate.) In Rrr, the ratio is enhanced at low T and 
high P, where quantum effects dominate, as we first expected in section ^ The ratio 
reaches 1.38 at T = 100 K and P = 1 GPa. Hence, the difference between the Yamomoto 
and the Miller approaches becomes rather large in the region of low T and high P. This 
tendency can be seen clearly in the contour plot of Rrr- On the contrary, in Rrp the ratio 
is reduced in the region where the quantum effect is strong. It is about 0.53 at T = 100 
K and P = 1 GPa. The contour plot shows the decreasing behavior of Rrp at low T and 
high P. 

Combining the (rr)- and (rp)-rate constants, one can calculate the ratio of the recom- 
bination rate constants, P^ec = kJ^JKL- The ratio is presented in Fig.||. (Note that the 
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ratio for the exchange process is given by Rexc = kJ^J^exc = Rrp-) The behavior of Rrec 
seems similar to Rrr and the ratio again reaches 1.38 at T = 100 K and P = 1 GPa. From 
the contour plot we can see that there is a small difference between Rrec and Rrr- 



3.2 Parabolic potential case 



The second example is a reaction which occurs under a harmonic oscillator potential. 
We suppose that the potential has a frequency u, the minimum point at Xq with its 
value Vq, and Sr and Sp are located symmetrically with respect to the minimum point 
(see Fig.|^). The propagator for a particle moving under the potential can be found by 
the path integral. The fiux-fiux autocorrelation function in the Miller approach is then 
calculated by (for details, see Appendix B) 



sinh(fi;/2) cosm 



(cosh(K/2) + cosm)^ sin^ u 



X exp 



(sinh^(/t/2) + sin^ m)3/2 (sinh2(K/2) + sin^ u)^/^ 

sinh(fi;/2)(cosh(/t/2) + cosu) 



(52) 



sinh^(K/2) + sin^ u 

where k = uj(3 and u = ut. Note that the correlation does not depend on the position of 
the minimum point explicitly. Similarly the correlation C*/ is obtained as 



Crr (t) 



sinh(K/2) cosM 



kCi 



(cosh(/«/2) — cosM)^sin^M 



X exp 



kCi sinh(K/2)(cosh(K/2) — cos u 



(53) 



sinh^(/t/2) + sin^ u 

Those correlation functions are shown by the dotted curves in Figs.|Ty and |TT| (we take 
Ci = 1.5 and n = 1.0 to illustrate the correlation functions clearly). In the limit tu, Vq — 
eqs.(^) and ( |55| ) are, of course, identical to eqs. (P^D and (P7|), respectively. 

After lengthy algebra, we can find the correlation functions in the Yamamoto approach 
(for details, see Appendix B). For example, for < m < tt, the correlations are expressed 
by 



Ke 



2\^nP'^ sinh kJ sin u(sinh^ k + sin^ u 



: exp 



kCi sinh ^(cosh n + cos u) 
2(sinh^ K + sin^ u) 



X 



Z_ cos ( —^Y. 



sin ( — 



(54) 
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and 



X 



2A/27r/?^ sinh sin ^(sinh^ k + sin^ u 



: exp 



kCi sinh /t(cosh n — cos m) 
2(sinh^ K + sin^ 



(55) 



where 



Z± 



sinh K + sin m ± sin 'u)(y sinh k + sin m ± cosh k sin -u 
1^1 ± cosm) sinh^ k — (cosh k — 1) sin^ u 



,1/2 



(56) 
(57) 



sin-u(sinh^ k + sin^ u) 

Note that in the hmit tu, Vq ^ the correlation functions approach those in the case of 
the free particle. Those correlation functions are also illustrated by the solid curves in 
Figs.0 and |n[ 

One can see from the figures that the correct behavior of the correlation function 
is quite complicated and it is divergent (like ~ 1/y/u) at u = 0,tt,2tt, ■ ■ ■. However, the 
Miller correlation function is smooth everywhere and it never diverges. All those functions 
are periodical because of the harmonic oscillator potential and, as expected, the integral 
of the correlation function over one period vanishes. 

The rate constants for recombination and exchange reactions are calculated by Laplace 
transforms of the correlation functions. In the Miller approach, the (rr)- and (rp)-rate 
constants are given by 



M 
rr 



-PVo 



AnP Jo 
X exp 



du 



sinh(K/2) cosM 
(sinh2(K/2) + sin2M)3/2 



(cosh(/t/2) — cosm)^ sin^ M 



(sinh2(K/2) + sin2M)5/2 



a kCi sinh(K/2)(cosh(/t/2) — cosu) 

^ ■ - 2/._ /o\ I „;„2 



K 



and 



M 
rp 



inP Jo 

a 



du 



sinh (k/2) + sin u 
sinh(/t/2) cosu 



(58) 



(cosh(/t/2) + cosm)^ sin^ u 



X exp 



u — kCi 



(sinh2(/c/2) + sin2M)3/2 (sinh2(K/2) + sin' u)5/2 

sinh(K/2)(cosh(K/2) + cosm) 



(59) 



K sinh'(K/2) + sin' u 

Similarly, we can obtain the Yamamoto rate constants. Because the expression of the rate 
constant is, however, lengthy, we do not write it explicitly here. (See Appendix B.) 
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Now we show our results of the parabohc potential case. In Figs.|T^ and |T^, the ratios 
Rrr and Rrp are illustrated. In the present calculation, we fix k to be 0.05, which means 
that the potential energy of the harmonic oscillator is much weaker (about 5%) than the 
typical thermal energy /3~^. We should note that the frequency u is varied so as to keep 
K = 0.05 at each T. If we set n to be smaller than 0.05 (for example, k = 0.01), the ratio, 
as it should, becomes close to that of the free particle case. 



In Fig.|12|, the ratio is again enhanced at low T and high P, which is similar to the 
result of the free particle case. The ratio at {T,P) = (lOOK, IGPa) is about 1.38. One 
distinct feature in the parabolic case is an enhancement of the ratio in the region of high 
T and low P. This can also be seen in the contour plot. In such a region, the power a 
appearing in Laplace transform for the rate constant is small, and hence the rate constant 
at high T and low P is more infiuenced by the correlation function at large t than that 
at other T and P, that is, the rate constant is considerably affected by the (second) 
complicated structure around m ~ vr in the correlation function (see Figs.|lO| and |ll]). This 
is the reason why the enhancement at high T and low P appears in the ratio. The ratio 
at {T,P) = (400K, O.lGPa) is about 1.11. In Fig.|TB| a similar tendency can be seen: the 
ratio is reduced at low T and high P {Rrp = 0.58 at {T,P) = (lOOK, IGPa)), which is 
similar to the free particle case, while it is enhanced at high T and low P {Rrp = 1.43 at 
(r,P) = (400K, O.lGPa)). 

Combining the (rr)- and (rp)-rate constants, we can calculate the rate constant for 
the recombination process; this is presented in Fig.jl^. In Rrec, the ratio is enhanced at 
low T and high P {Rrec = 1-38 at (T, P) = (lOOK, IGPa)), while it is reduced at high T 
and low P {Rrec = 0.38 at (T, P) = (400K, O.lGPa)). This behavior is also seen in the 
contour plot. 
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4 Summary and Conclusion 

The exact quantum mechanical expression for thermal reaction rates can be formulated 
by the linear response theory,!* I which Yamamoto first discussed in the early 60's.i Later, 
in the early 70's, Miller et al.S i have independently proposed a more convenient way to 
perform numerical computation, which can provide the exact rate constant in the limit 
that the dynamics of the system is extended to t ^ oo. 

We have studied the difference between the two approaches in thermal reactions which 
involve exchange and recombination processes. Because the rate constants in those reac- 
tions are calculated by Laplace transforms of the flux-flux autocorrelation functions, the 
results evaluated by the two approaches are different. In this paper, we have considered 
two solvable cases, i.e., the free particle and parabolic potential models, to demonstrate 
the difference intuitively. We have found that the shapes of the correlation functions 
are quite different in the two approaches and that the difference of the rate constants 
appears in the region where quantum effects dominate. In both the free and parabolic 
cases, the rate constant for recombination in the Yamamoto approach is larger than that 
in the Miller approach; the enhancement becomes about 40% at low temperature and 
high pressure. 

In conclusion, the Miller method is certainly an economical and powerful tool to 
perform numerical calculations for thermal rates of realistic reactions. However, it may 
underestimate the rate constants for recombination and exchange processes in the region 
where quantum effects dominate, because of neglecting the A integral appearing in the 
Kubo identity. 
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Appendix A — free particle case 



The matrix element of the flux operator Fi{i = r, or p) in the coordinate representation 
is given bj 



{u\Fi\u') = -—[S'{u - s,)6{u' - s,) - 6iu - s,)6'iu' - s,)]. 



2im 



(60) 



Using this expression, one can evaluate the flux-flux autocorrelation function in the Miller 
approach as 



C^(T) = tr(F,(0)e^*^^Fp(0)e 



it*H ; 



-itcH\ 



2m2 



3? 



d_ 

du 



{u\e''"'^\u'Y^iu\e-''''-\u') 



du' 



ouou' 



(61) 



J U=Sr.U '■ 



where tc = t — ij3/2 and 3? stands for taking the real part. On the other hand, the 
Yamamoto correlation function is given by 

Cr,{T) = ^tr(e-^^[M^.),F,(t)]), 



1 f°° 
— I du 

mp 



i{u\e-''''^\u'y^{u\e-''''\u') 
ou 



- i(u\e-'^'\u')^{u\e~'^'^\u'y 
du' 



(62) 



where tjs = t — ijS and stands for the imaginary part. 

With use of the path integral technique,^ the propagator for the free particle at finite 
f3 is calculated by 



{ule-'^^'^lu') 



2m (t2 + /?2)iA 



exp 



where 



cos^^ 



== and sin 9 



m{it - (5) 2 

/5 



(63) 



(64) 



We can easily calculate the correlation functions using those expressions. The final 
results are explicitly presented in eqs.(|36D-(^0D. 
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Appendix B — parabolic potential case 



The propagator for a particle moving under the harmonic oscillator potential (at finite 
/?) which has the minimum point at (a;o, Vq) (see Fig.|^) is evaluated byEl 



{u\e 



-iHta 



V V27ri 
X exp 



f mLj\ cosh Ksinu + isinh Kcosu _^y^ 



tmuj 



sinh K, + sin u 
f sin u cos u + i sinh k cosh k ' 
\ 2(sinh^ K + sin^ u) 



/ cosh K sin m + z sinh k cos u 



imuj 



uu 



iVot 



(65) 



\ sinh K + sin u 
where k = ujS, u = ut and tp = t — i(3. 

Because it is easy to calculate the Miller correlation function using eqs.(0) and (|65D , 
we do not show its derivation here. The final result is given in eqs.(^) and (|53|) . Instead, 
we explicitly present the Yamamoto correlation function. It is a periodical function and 
the period is An (see Figs.|10|and |1T]). Then, we divide it into 4 parts. Using eqs.(^) and 
(|65|), we find: 

(1) for < M < vr, the results are given by eqs.(0) and f|5l 

(2) for TT < M < 27r, 



2v^vr/52 sinh kJ \ sin u| (sinh^ k + sin^ 



exp 



u) 



kCi sinh ^(cosh k, — cos u) 



2(sinh K + sin u) 



X 



Za. cos ( — 



and 



Z_ sin 



fKC_ 



V 2 



(66) 



2-\/2vr/5^ sinh kJ \ sin -ul (sinh^ k + sin^ u] 



exp 



kCi sinh ^(cosh n + cos u) 



2 (sinh K + sin u) 



Z_sin('^n 



Z I cos f -— , , , 

. V 2 +y V 2 

where Z± and Y± are defined by eqs. (|56D and 

(3) for 271 <u< Svr, 



(67) 



C„{t) 



2v^vr/52 sinh kJ sin u(sinh^ k + sin^ u] 



exp 



kCi sinh ^(cosh k — cos m) 
2(sinh^ K + sin^ m) 



X 



.Z_cos('^F_ 



(6^ 
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and 



cut) 



ne 



-/3Vb 



2^/2n(3'^ sinh kJ sin u(sinh^ k + sin^ u 



: exp 



-Z_ cos ( — — !-y4 



Z4. sin ( — — !-y4 



(4) for Svr < M < 47r 



ne 



X 



2-\/2vr/?^ sinh /ty | sin (sinh^ k + sin^ 

cos sin 



exp 



and 



KC 



2y/2Tt(3'^ sinh | sinM|(sinh^ k + sin^ -u' 



exp 



-Z< cos ( — ^F-i 



Z_sin('^n 



kCi sinh fi;(cosh k + cos u) 
2(sinh^ K + sin^ u) 



(69) 



kCi sinh ^(cosh k — cos m) 



2 (sinh K + sin 



(70) 



kCi sinh ^(cosh k + cos u) 
2(sinh^ K + sin^ u) 



2 --7 ---V 2 --yj 

As shown in Figs.|l^ and 0, there is a symmetry property: C^(0 < u < vr) = 
-C^^(27r < M < Stt) and C^^(7r < m < 27r) = -C,^(37r < u < 47r) (i = r or p). Thus, 
it is enough to calculate the correlation function in the region of < m < 27r. The rate 
constants are calculated by Laplace transforms of those periodical correlation functions. 
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Figure captions 



Fig. 1: Sketch of a potential surface in one- dimensional reaction versus the reaction 
coordinate q. 

Fig. 2: One-dimensional schematic diagram of the interaction potential for the A + BC 
AB + C reaction. The compound region, ABC, is bounded by the dividing surfaces on 
reactant (s^) and product (sp) sides. 
Fig. 3: Same as Fig.^ but for the free particle case. 

Fig. 4: Correlation function C„(t) for the free particle. The dotted curve is for C*/(t), 
while the sohd curve is for C^(t). We take Ci = md? /2[3 = 1.0. 

Fig. 5: Correlation function Crpit) for the free particle. The dotted curve is for C^p(t), 
while the sohd curve is for C^pit). We take Ci = md'^/2(3 = 1.0. 

Fig. 6: Ratio of the (rr)-rate constants in the free particle case (top) and the contour 
plot (bottom). In the contour plot, the top dotted curve corresponds to Rrr = 1.05 and 
the other curves are plotted at intervals of 0.05. The bottom and right dotted curve is 
thus for R„. = 1.35. 

Fig. 7: Same as Fig.^ but for the (rp)-rate constant. In the contour plot, the top dotted 
curve corresponds to Rrp = 0.95 and the other curves are plotted at intervals of 0.05. The 
bottom and right curve is for R^p = 0.55. 

Fig. 8: Same as Fig.^ but for the rate constant for the recombination reaction. In the 
contour plot, the top dotted curve corresponds to Rrec = 1-05 and the other curves are 
plotted at intervals of 0.05. The bottom and right curve is for R^ec = 1.35. 
Fig. 9: Same as Fig.^but for the parabolic potential case. The minimum point is located 
at (xo, Vo). 

Fig. 10: Correlation function Crr{t) for the harmonic oscillator case. The dotted curve is 
for C^^(t), while the solid curve is for C^(t). We take Ci = 1.5 and k = 1.0. 
Fig. 11: Correlation function Crp{t) for the harmonic oscillator case. The dotted curve is 
for C^p(t), while the solid curve is for CZ,(t). We take Ci = 1.5 and k = 1.0. 
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Fig. 12: Ratio of the (rr)-rate constants under the harmonic oscillator potential with 
K = 0.05 (top) and the contour plot (bottom). In the contour plot, the dotted curve 
connecting (T, P) = (lOOK, 0.14GPa) and (380K, IGPa) corresponds to i?„ = 1.05. The 
top, left dotted curve connecting (T, P) = (340K, O.lGPa) and (400K, 0.14GPa) is also 
for R„. = 1.05. The other curves are plotted at intervals of 0.05. The bottom and right 
dotted curve is thus for P^r = 1-35. 

Fig. 13: Same as Fig.|l2| but for the (rp)-rate constant. In the contour plot, the dotted 
curve connecting (T, P) = (lOOK, 0.33GPa) and (210K, IGPa) corresponds to Rrp = 0.95, 
while the dot-dashed curve connecting (T, P) = (270K, O.lGPa) and (400K, 0.19GPa) is 
for Rrp = 1.05. The other curves are plotted at intervals of 0.1. The bottom and right 
curve is for Rrp = 0.65. 



Fig. 14: Same as Fig.|12| but for the rate constant for the recombination reaction. In 
the contour plot, the dot-dashed curve connecting (T,P) = (120K, O.lGPa) and (400K, 
0.93GPa) corresponds to Rrec = 1-05, while the dotted curve connecting (T,P) = (275K, 
O.lGPa) and (400K, O.lSGPa) is for Rrec = 0.95. The other curves are plotted at intervals 
of 0.1. The bottom and right curve is for Rrec = 1-35. 
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Figure 13: K. Saito 
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